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TURBULENT BOUNDARY LAYERS WITH SECONDARY FLOW 


E. Grushwitz 

The previous publications on friction layers or boundary 
layers^- concern almost exclusively bidimensional flow processes, 
but naturally the processes with rotation symmetry have also been 
studied. This study represents an attempt to learn experimentally 
something about the processes in a boundary layer with three- 
dimensional flow, and specifically the boundary layer is studied 
on a plane wall, along which a flow occurs, whose potential flow 
lines are curved in planes parallel to the wall. A pressure 
gradient directed perpendicularly to the flow line on its concave 
side and lying in the plane of curvature corresponds to a flow 
line curvature. Since the static pressure in the boundary layer 
is mainly the same as in the limiting potential flow, the flow 
line in the boundary layer, because of the lower velocities occur- 
ring there may be curved more strongly than the potential flow lines 
or assume a different direction from them. The flow lines in the 
boundary layer are consequently deflected more toward the inside, 
the closer they are to the wall. The name secondary flow has been 
given to the component arriving there perpendicular to the direc- 
tion of the potential velocity. 

In view of the difficulties arising in the comparison of 
spatial flows according to their mathematical treatment, and since 
very little is known up to now about the turbulence mechanism. 


iFor the region adjacent to a body surrounded by flow, in which the 
decrease of velocity to zero takes place, the author uses, to satis- 
fy the wish of Professor Prandtl, the word "friction layer" instead 
of the "boundary layer" commonly used in the literature. The boun- 
dary layer will designate the laminar component of the friction 
layer. 


- *Numbers in the margin indicate pagination in the foreign text. 
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there is hardly any prospect of mastering even approximately these 
complicated processes under the theoretical aspect within a 
foreseeable time. The measurements described in this article are 
therefore carried out to obtain an indication of the order of 
magnitude of the effects occurring in this connection. 

I have to thank Professor Prandtl for suggesting this work and 
the Cooperative Aid Association of German Science for making it 
possible. 

1 . Experiments 

1. The experimental channel . The experimental equipment 
(see Fig. 1) consisted of a channel of rectangular cross section, 
which had first a straight piece 1 m long. This was followed by 
a channel piece curved by 90 the diameter of curvature of the 
internal channel wall was 45 cm, that on the external wall 90 cm. 
The boundary layer of the upper plane channel was studied. For 
this purpose in this channel wall, consisting of 5 mm thick iron 
plate, holes were provided, whose position is apparent from Figs. 

1 and 5 (circles in Fig. 1). These holes y/ere ' provided with 
thread, in which it was possible to screw a shifting device called 
the measurement equipment used to study the boundary layer. The 
unused holes were sealed with screw plugs, which closed flush 
with the internal channel wall. In the places (points) shown in 
Fig. 1, holes were also provided to measure static pressure. Air 
was blown through the channel by means of a blower in the direction 
indicated by the arrow. Between the blower and the experimental 
section a rectifier, several sieves and a 2.5 m long stabilizing 
channel were installed. The pressure in the stabilizing channel 
and therefore the wind velocity were kept unaltered through a 
Schrenk regulator^. The blower was the same as used earlier by 
the author to measure boundary layers in a plane flow^. Just as in 


2 

Compare 0. Schrenk, Ing-Arch. 1 , 350 , (1930) . 
^E. Gruschwitz, Inq.-Arch 2 , 321 (1931) . 
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the case of those measurements, here too a slit was provided in 
front of the measurement plate, through which the friction layer 
coming from the wall of the stabilizing channel was able to flow 
out of the channel, so that a new boundary layer began on the 
measurement plate. As shown by Fig. 1, the experimental channel 
had in its front part a small extension and from then on its 
cross section remained the same (45 cm. wide,. 50 cm high). The 
extension was installed for the purpose of allowing the boundary 
layer at the measurement plate to become turbulent as soon as 
possible. 



Figure 1. The experimental channel and 
measurement plate (measurement plate seen 
from below) . 

2. The measurement instrument . In the boundary layer to be 
studied there are mainly two velocity components: specifically 

the two components parallel to the wall. The third, perpendicular 
to the wall is so small as compared with the other two, that it 
may be considered negligible for the measurement. Therefore 
measurement equipment had to be used, which allowed the measurement 
of these two components or with which the value of the velocity 
of the wind and its direction could be measured. As such an 
instrument a thin cylindrical tube was used, which could be 
shifted perpendicularly to the wall in the flow, through a hole 
in the wall on which the boundary layer was to be measured. In 
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this tube three holes were provided near each other in a plane per- 
pendicular to the tube axis. The tube was turned until both 
external holes showed the same pressure, by which the direction 
of the velocity of the wind was determined. Then the total pres- 
sure was measured with a third hole. The free end of the cylinder 
tube was closed with a hemisphere, whose diameter was equal to 
the cylinder diameter. 

To obtain a survey of the variation of 
pressure at the end of the tube, a similar 
larger one was first built, whose diameter 
was 11.8 mm. In it pressure measurement 
holes were provided in the places indicated 
in Fig. 2. The cylindrical tube was brought 
into an air flow in such a way that the 
generator on which the pressure measurement 
holes lay became the straight line of stag- 
nation. In Fig. 2 the measured pressure q z , 
refer to the pressure head q Q of the undis- 
turbed flow as plotted over the distance 
from the tube end. The blowing speed was 
in this case 20.5 m/sec. 

The position of the holes in the tube 
used finally for the measurements may be 
seen in Fig. 3. The pressure at both external holes was sent up- 
wards through the measurement tube through thin steel tubes (internal 
diameter 0 . 4 mm, external diameter 0 . 6 mm) through the measurement 
tube, while the central hole ended directly in the measurement 
tube. The external holes were connected with an oblique U-shaped 
tube, the central one to a Prandtl manometer. Both thin steel 
tubes gave rise unfortunately to fairly long adjustment time. 

The central hole did not show fully the total pressure head, 
therefore a calibration factor must be taken into consideration. 

If the indication of the cylindrical tube is designated by q z and 



Figure 2 . Pressure 
distribution along 
the straight line 
of stagnation of a 
cylindrical tube 
with hemispheric 
end with flow 
perpendicular to 
its axis. 
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the pressure head by q, then, as shown by compari- 
son with a pitot tube, q=1.028q z . 

If the total pressure is measured in a boundary 
layer with the described cylindrical tube, it may 
be expected, that the disturbance caused by the 
measurement tube in the boundary layer, is greater 
than when the measurement is carried out in a 
pitot tube, since in the latter only the mouth of 
the tube is located at the point to be studied. 

In particular in the cylindrical tube a distur- 
bance will be noted by the fact that in the 
boundary layer the pressure along the stagnation 
line of the tube is not constant, but becomes smaller near the 
wall as a result of the decrease of the velocity. This circumstance 
will presumably cause in front of the tube a flow toward the wall, 
so that the flow lines of the boundary layer flow are deflected 
somewhat toward the wall, which led to the fact that the velocities 
in the boundary layer were measured somewhat too large, especially 
wherever the velocity gradient is large, that is, near the wall. \ 
To establish the differences between the indication of the cylindri- 
cal tube and that of a pitot tube, a velocity profile of a plane 
boundary layer was measured consecutively with both instruments. 

The graph in Fig. 4 shows that the cylindrical tube actually gives 
a somewhat too large velocity near the wall, on the other hand the 
velocity is somewhat too low in the external regions of the boundary 
layer. As may be seen from the graph, these errors are however not 
very large and they had to be accepted, sinde it hardly seems 
possible to introduce a correction, since its order of magnitude 
depends obviously from the corresponding variation of the velocity 
profile. 

A measurement tube was attached to a shifting device, which 
makes it possible to shift it along the direction of its axis and 
to turn it around its axis. It was possible to measure the shift 
with a precision of 0.01 mm, and the rotation with a precision of 



Figure 3 . The 

measurement 

tube. 
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10'. This shifting device was screwed 
into the above mentioned, hole located 
in the channel wall. 

3. Measurement results. Figure /358 
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Figure 4. Comparative 
measurement between pitot 
tube and cylindrical tube. 


Key: 


a. cylindrical tube 

b. pitot tube 


5 shows the measurement plate looking 
from below. The points at which the 
velocity profiles were measured are in 
the intersections of the lines indicated 
Roman and Arabic numbers. The lines 
designated with Roman numbers consist 
each of a straight line section and a 
curved arc, those designated by Arabic 
numbers are straight lines orthogonal 
to them. In Fig. 5 some flow lines on the potential flow are shown, 
plotted according to the measurements of angles with the cylindri- 
cal tube. In the shaded region there was no longer any potential 
flow; the boundary layers of the lower and upper channel wall ran 
into each other there. In the last portion of the channel only a 
small region with potential flow occurred. On the internal wall 
of the channel curvature the flow was torn off and specifically 
only in the central portion of the curved channel wall, whereas it 
was applied in the corners. This behavior which may appear re- 
markable at first glance is obviously to be attributed to the 
effect of the secondary flow, which continues from the upper and 
lower channel wall to the internal wall and causes a combined flow 
of the higher delayed air layers in the center of this wall. 


The variation of the static pressure of the measurement plate 
is shown on Fig. 6. There the static pressure p refers to the 
total pressure g Q of the potential flow and measured on the pres- 
sure holes against the external chamber, is plotted over a series 
of cross sections. 


Let the velocity of the potential flow be designated by U, 
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Figure 5. Measurement 
plate seen from below. 

The measurement points 
are at the intersections 
of the curves or straight 
lines designated by 
Roman or Arabic numbers. 
The solid line curves 
are flow lines of the 
potential flow. 


Figure 6. The static 
pressure on the measure 
ment plate. 


the component of the boundary layer flow in the direction of the 
potential flow by u, the one perpendicular — by v. In so far as 
we are dealing with turbulent flows, we always refer to averages 
in time. Figure 7 shows as an example the velocity profile measured 
along the central line of a channel (line III in Fig. 5) , and 
specifically the velocities u and v are referred to the potential 
velocity U occurring at the measurement point concerned, and 
plotted over the distance z from the wall. Figure 8 shows the 
variation of a velocity at some point of line III as a polar dia- 
gram, while the heights z are indicated in millimeters at individual 
points. 
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Figure 7. Velocity profile along the central line of 
a channel (Curve III in Fig. 5) . The numbers on the 
profiles correspond to the designations in Fig. 5. 



Figure 8. Polar diagram of the velocity at 
the measurement points 6, 9, 12, 15, 18 and 
20 of Line III. The distances z from the wall 
are indicated in millimeters at the individual 
points . 


Let Uq designate the potential velocity, corresponding to the 
static pressure p=0. Therefore we have JgaU^go- The value U Q is 
essentially the velocity with which the potential flow leaves the 
channel. In the described measurements we have U 0 =13.2 m/sec. 
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II. Application of the Pulse Equation to the Boundary Layer with 
Secondary Flow 

1. Transformation of equations of motion . The hydrodynamic 
equations of motion will be converted. as follows: let us imagine 

a flow line of the potential flow projected perpendicularly to 
the wall on which the boundary layer is located. Let x be the arc 
length along this projection, let y be the distance of the normal 
to it lying in the wall, and z the perpendicular distance from the 
wall. The equations of motion converted in this system apply in 
the vicinity of a cylindrical surface lying perpendicular to the 
wall and placed through the potential flow line concerned. Let 
X, Y, Z be the coordinates in a rectangular system of ordinates, 
whose (X, Y) plane coincides with the wall. It is apparent from 
Fig. 9 that the following equations of conversion apply: 

x 

X = J cesa dx — y sina, 
o 

Y = Fi-f ycosa, 

Z = s. 


(1) 

l 


Here Y^ is the ordinate of the projection of the potential flow line 
and a the angle between the tangent to this curve and the X axis. 
Both and a are here to be considered as functions of the arc 

dec 

x. The curvature of the curve is designated by k. We have k= 


By differentiation we obtain from (1) 


(I— yk) cosot, 

\ y = (* — y k ) sina > 

dZ 

dx ~ °’ 


ex 

dy 
dY 
d y 
dZ 
dy 


= — sin a , 


= cos a, 


= o, 


dx 


dz 

dY 


dz 

dz 

dz 


o, 


= I 


We obtain thus for the, functional determinant D 
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= I — yk 



F.igure 9. 
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and consequently for the derivatives of the 
inverse functions: 
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( 2 ) 


The equations of motion in the rectangular system of coordinates 
X, Y, Z are: 


e[ u °TX + v °lT + “'»4ar) 
e(« 0 |5 ? + «( 


8a, » . 8tt '<A 
0 ay ^ u ° az j 


box , 

a txy 

4_ 

brxz 
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ay 
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ay 1 

az 

i 

ax 

ay 

a az , 

, Staz 

4- 

aTyz 

bp 
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i 

ay 
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Moreover we have the equation of continuity: 

dtt 0 , dv 0 , _0w_o_ n 
ax + 8Y ‘ r 0Z 

Here the velocities in the X, Y, Z directions are designated by 
u D , v Q and w Q , the density by a, the static pressure by p and the 
matrix of the (apparent) friction stresses by 

( ox x x v t xz 
Txy oy r yz 
T xz 1YZ OZ 



We obtain with formula (2) after some calculation, the transformed 
equations of motion: 


i 


I 
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( u dn . du . du\ k 

^ \ i — yk d x V dy $ z ) Q UV ~ J — yk ~~ 


C* Ox 


yk 

&?xy , d T x z 2 T x y k 

' d z I — yk 


d x i — y k ' dy 
( u dv . dv , dv \ x « h 

Q \T-yk T7 +v Ty + w ~dl) + e u 7 -Jk = 


_l dp 
i — y k 0 x ' 


__ da y 


dT; 


xy 


dy i — y k d x 


tijyz 

dz 




((Tx (jy)k 


( dw , duo dw\ dox , I drxt , 

e[ u J7 + V Ty + w u) = -8- + T-y-k-Fr + 


i — yk 


d y 


dp 
dy 1 

T yxk 
l yk 


d p 
dz 1 


and moreover the equation of continuity 

1 1^1 dlV V k n 

i — ykdx'dy-^dz i — y k ~~ 


Here u, v, w are the velocities in the s, y, and z direction and 

j Ox *xy T x z 
i Txy Oy Tyz 
\t xz T yl Oz 



the matrix of the friction stresses in the (x, y, z) system of 
coordinates. The transition from the friction stresses in the 
(X, Y, Z) system to those in the (x, y, z) system is here given 
by the matrix equation 


( ax txy t X z 
txy oy tyz 
tx z tyz oz 


X 


C cos(arA), cos(yA), cos (oA)\ / o x T xy t xz \ 
cos (xY), cos(yF), cos(cF) M T xy o y T y , I x 
cos {xZ) , cos (yZ) , cos ( zZ)j \ r„ T yx a 2 ) 

( cos(*A r ) ( cos(*F), cos(xZ)\ 
cos(yA), cos(yF), cos (xZ) J = 
cos (sX) , cos (s Y) , cos (zZ) ) 


'cos a, —-sina, o 
sina, cosa, O 
o Ol, 


/ O x 

Txy 

Txz \ 

t cosa, 

sina, 

°\ 

Txy 

Oy 

Tyz 

I — sina, 

cosa, 

0 

w 

Tyz 


V ° 

0 

1/ 


where (xX) designates the angle between the x and X direction, (yX) 
the angle between the y and X direction etc. 


2. The pulse equation . We will now integrate the first of 
equations (3) for y=0 over 2, assuming that the static pressure 
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is independent of z inside the boundary layer. Of the friction 
terms here only 3t xz /9z can be maintained, since in the boundary 
layer the other friction terms are small as compared with this. 
With this omission and for y=0, the first equation of motion is: 


»£ + *7f + "Tr-*'' k = 


du 


Bt X z T. dp 

~dz q dx 


( 4 ) 


and the equation of continuity 


du 

~dx 


. dv , dw y 

+ jy +sT~ vk 



(5) 


Equation (4) is converted into 


1 
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d u z , duv 
d x * cy 


dv . d u w 

U u b 

d y 1 


dz 


d w 7 

- u 0 uvk 

dz 


I 

(j dz 


i dp f 

u & x * j 


In this equation we replace -r— by the expression which is obtained 

d Z 

for it (5) and we now have 


f 



du v 

~dY 


+ 


duw 

dz 


2 U V k 


_I_ dixz I dp * 

Q dz Q dx ’ j 


This equation integrated over z gives 

c e c 

■ -^-fuUz + 4-fuvdz+UlV—2kfuvd2 = — -r;l— c -^t. 
8 x J 1 8 y J 1 • J Q * 0 8 x 

o. o ...o 


7364 


( 6 ) 


Here c is a constant length, which is so large that for z=c, 
potential flow is present practically everywhere (therefore t X z=0) / 
i s the value of x X z (shear stress in the x direction) on the 
wall, and U and W designate the x and z components respectively of 
the potential velocity. We obtain for W by integration of equation 
(5) over z 

c e . c 

W = kjvd,—^fud^jvi=. (7 

0 _ 0 0 

According to the Bernoulli equation we have 
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The Roman and Arabic numbers, designating the measurement points, 
correspond to the designations in Fig. 5. Xj is the distance of 
the measurement points from the front edge of the measurement plate, 
measured on the lines designated in Fig. 5 with Roman numerals. 
Distances between the lines are 7 cm. 4> 0 is the angle which is 
formed by the boundary flow line on the wall with the direction of 
the potential velocity. The angle 4> 0 is obtained by extrapolation 
of the angle measured with the cylindrical tube. 


14 



REPRODUCED FROM BEST 
AVAILABLE COPY 


‘ ~ i' 


Kc ray also write for it 


‘ */ V ' /'/*./; — r //* 4 

7 t 9 / n; 


Moreover we obtain for the curvature of tho potential flow line 


as obtained from the second of equations (3) together with the 
Bernoulli equation. 

If (7), (0) and (9) are introduced into (G) , we obtain finally 

* * i * 

/;/'•-» »J:~- * ' I r~ » ./:4. * ■ f-f-u-i-Ji— ?. /./ -« i</: (10) 

• * / < ( / | / / i v ' M 

« *# ‘ « 

Mow the following designations rvay be introduced: 

; - l ‘ . 


Kith these designations we obtain frora (10) 


‘ ■*« • • . . 

"7* /> • l'i or 

« « * II 1 


< r «. / * 


This equation nay bo called the "pulse equation". It is a generalization 
of the eejuation valid for friction layers in a plane flow, which 
was first derived by Karr.an by applying the pulse flow to the 
friction layer or boundary layer 1 *. 


For the sake of completeness we should also indicate the 
^Th. v. Karmen, Z. ar.aew. Hath. Mgch. 1, 233 (1521). 
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corresponding integration valid for the y direction, although it 
is not used hereafter. It may be derived .in the same, way as 
equation (11) from the second part of Eq. (3) for y = 0 and dis- 
regarding the friction term except for ■ — , and becomes i 

o Z 


d &y , d (<3 y &x y) 

d y ' 


dx 


U dy 


BU la n C X . 2 dU 

(& y — &x — <5*) + jf TfJ (°y ' 


T (i0 

y* 

nU 2 


Here the designations 


&y = -fj-Jvds 
0 w 0 


are further introduced and is the value of Ty Z of the wall. 

3. Application of the pulse equation to the measured velocity 
profile . In the above Table the values of 0 X » 0y/ ©xy* ^x an ^ ^y 
obtained from the measured velocity profile are indicated. To give 
a survey of the order of magnitude of the individual elements of 
equation (11) the latter are also given in the Table. The deriva-^ 
tives in the three terms on the left side of equation (11) were 
established graphically from the values obtained with the measure- 
ments. It was possible to determine them only with rather high 
uncertainty, in particular, 30 X y/3y was extremely hard to establish. 
Wherever values are not indicated in the Table, it was impos- 
sible to determine the derivatives. The values obtained by addi- 
tion of the three terms for /oU 2 are to be considered only as 

estimated values. It is apparent from the Table that the term 
representing the effect of a secondary flow in (11) is of the same 
order of magnitude as the others. 


For comparison the Table also gives the wall shear stresses, 
which would be present in a boundary layer in plane flow of the 
same thickness (same 0 X ) for the same potential velocity. For the 
Reynold numbers occurring here, the assumption of a velocity 
distribution in which the velocity is proportional to the 7th root 


/365 
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of the distance from the wall, provides a sufficient approximation 
for these shear stresses. With this assumption we obtain for the 
shear stress t 0 of the wall^ 

Jv * =0,01255 

It is apparent from the Table that the. values of shear stress 
obtained from the best measurements are spread around those in the 
plane flow. 

Summary 

The article describes measurements with an experimental equipment used to study 
the boundary layer on 'a plane wall, along which the flow occurs, whose potential 
flow lines are curved in plane parallel to the wall. According to the equation 
frequently applied to boundary layers in a plane flow, wich is usually obtained 
there by using the pulse law, a generalization is derived which is valid for 

boundary layers with spatial flow. The wall shear stresses were calculated 
with this equation. 


^Compare v. Kdrmdn, Loc. cit and E Gruschwitz, Inq.-Arch. 2 , 341 
(1931). 
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